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ABSTRACT 

We discuss a special case of formation of axisymmetric shocks in the accretion flow 
of ideal gas onto a Schwarzschild black hole: when the total energy of the flow is 
negative. The result of our analysis enlarges the parameter space for which these 
steady shocks are exhibited in the accretion of gas rotating around relativistic stellar 
objects. Since keplerian disks have negative total energy, we guess that, in this energy 
range, the production of the shock phenomenon might be easier than in the case of 
positive energy. So our outcome reinforces the view that sub-keplerian flows of matter 
may significantly affect the physics of the high energy radiation emission from black 
hole candidates. We give a simple procedure to obtain analytically the position of 
the shocks. The comparison of the analytical results with the data of ID and 2D 
axisymmetric numerical simulations confirms that the shocks form and are stable. 
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1 INTRODUCTION 

The presence of shocks in accretion flows is a relevant fact 
in the theory of accretion discs. Shocks are a very effi- 
cient mechanism to produce high energy emission of radi- 
ation and to induce strong acceleration of particles. Spiral 
shocks in accretion disks in binary systems h ave been stud - 
ied since long time, theoretically by Spruit (|SErui3[l983), 
Chakrabarti llChakrabartil Jl 990f) and with both 2D and 
3D numerical simulations (ISawada et al ]|l987l : iTaam et alJ 
Il99ll: iBlondin et al.ll995l: lYukawa et alJll997ll . Furthermore 
there is increasing observational evidence of their significant 
role in determining the accretion of gas in binary systems 
iHachisu et al.ll2004) . However the strength of this kind of 
shocks is usually weak. 

A different scenario for sho ck producti o n ha s 
been proposed by Chakrabarti dChakrabartil ll99(J) . 
It ha s been widely tested with numerical s i mula- 
tions (ICh akraba rti fc MoUenil Il993l: iMoheni etHI Il994t 
IRvu et 5 . 1996). In this view, an axisymmetric flow of gas, 
rotating with an angular velocity less than keplerian (we 
will call it "sub-keplerian" flow), and with small viscosity, 
can produce a steady shock. This fact is mainly due to the 
deceleration of the radial fall and consequent heating of the 
gas close to the centrifugal barrier. Although the origin of 
such flows seems unclear there are arguments to support the 
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hypothesis that they really exist. In galactic nuclei the am- 
bient gas feeding the central black hole may be quite well 
symmetric and with sma ll rotation as shown by observations 
JCohen fc Rvzhovll99i1) . On the theoretical side, taking into 
account also the role of the viscosity, a detailed study of the 
requireme nts for their origin from a keplerian disk has been 
produced iChakrabartilll99d) . These flows may explain very 
well the two frequencie s in the QPO phenomeno n in galactic 
black hole candidates dChakrabarti et al]|2004l) . So even if 
their direct proof of existence is still pending, we think it is 
useful to study their properties and characteristics. 

In any case it seems to us that the phenomenon, de- 
spite of its basic physical straightforwardness, has not been 
adequately considered. Maybe some "unconscious" ques- 
tions are raised against their effective existence: 1) These 
shocks appear when the flow has well defined parameters, 
derived from the analytical treatment. 2) In the case of non 
ideal gas (i.e. if a real cooling, an heating or conduction 
are included) the analytical treatment is very difficult. 3) 
Even in the simple ideal gas case, the standard procedure 
JChakrabartill99Cl) to derive a solution is rather compli- 
cated: one has to compute the sonic point location and then 
compute the boundary conditions on derivatives of the ra- 
dial speed, sound speed etc. at the sonic point itself, leading 
to cumbersome formulae. 4) At the end of the whole process 
the analytical solution is numerical and specific for the set 
of the chosen parameters. So no simple general formula, like 
the Shakura- Sunyaev disk structure formulae, is available. 
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Another doubt is that the parameter space seems to be 
Hmited or too compelhng. Roughly speaking, people might 
say that the sources of gas do not know the boundary condi- 
tions at the sonic point or at infinity, they spread out their 
gas as they like. So, for example, in bound systems, the to- 
tal energy will be negative and, up to now, only positive 
energy solutions have been discussed and simulated . Even 
standard ADAF models (see lBlandford fc Beeelmanl lll999l) 
and references therein), that are similar to the mentioned 
sub-keplerian flows, require positive energies to flt observa- 
tional data. So what about negative energy flows? Here, we 
explain that, even in this case, there are solutions admitting 
permanent shocks. Obviously this does not mean that in 
every sub-keplerian flow there is a standing shock at work. 
But when the physical parameters of the accretion (essen- 
tially the rotation amount and the thermal energy content) 
are appropriate, then the shock process occurs and produces 
the relevant phenomena described in previous references. 

The arguments of the paper are given as follows: in 
section 2 we revise the general accretion formulae for ideal 
gas with negative energy, in Section 3 we present the results 
of numerical simulations and in Section 4 we make our final 
discussion. 



is given by: 

^ + + + ^ + =E (3) 

We will look for isentropic solutions and then we will see 
how to connect solutions with different entropies by shocks. 
For a flow with constant entropy S and with specific heat at 
constant volume Cv we may exploit the adiabatic relation: 



p = e"c.(T-i) {^\ (4) 

that allows us to eliminate the density p from the equa- 
tions and to use the sound speed. 

So we have two unknown quantities Vr , a and two equa- 
tions. Resolving for a from the Bernoulli relation, using the 
radial Mach number M — — — , and putting all terms into 
the continuity equation, we have the following implicit solu- 
tion for the Mach number: 

M = -47ryl(r)/(M) (5) 

where f{M) is the Mach function which depends only 
on the M values, given by: 



2 ACCRETION OF IDEAL GAS WITH 
NEGATIVE ENERGY 

To explain the shock formation, let us resume the very ba- 
sic physical ingredients of the phenomenon. We assume that 
an axisymmetric fiow of an inviscid gas is falling from very 
large distance onto a black hole and has reached a steady 
state regime. We will adop t the Paczyhski & Wiita potential 
JPaczvhski fc Wiital[l98ol) to mimic the general relativistic 
physics. The basic physical effect can be easily understood 
in terms of "classical" physics and it is well known that the 
Paczyiiski & Wiita force reproduces many relativistic effects 
with high accuracy. As we said in the introduction, to obtain 
the steady state solution we may integrate the differential 
equations for mass, momentum and energy conservations. 
In this case one has to start the space integration from th e 
sonic point as explained in the work bv IChakrabartil il99d) . 
However in the case of inviscid flow it is easy to find an alge- 
braic implicit solution, since (in this case) the total energy 
is conserved. The Bernoulli theorem is valid and it can be 
exploited to close the system of equations and to find the 
solutions. For the 2D case (with real Z extension of the disk) 
we assume that the gas falls down in a condition of vertical 
equilibrium, with essentially zero vertical speed. 
The mass conservation equation is given by: 



M 



-ATvrH pvr ~ const 



(1) 



where H is the half thickness of the disk. 
The vertical equilibrium hypothesis gives the well known 
expression for the half disk thickness: 



H = 



(2) 



We could use a more accurate formula, but the results 
are essentially unchanged while the formulae are more com- 
plicated. 

According to the Bernoulli theorem the energy equation 



/(M) = q- 



M 
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where g = (7 — 1)'''""^^ 

This function (we name it the "Mach" function) has a 
maximum at Mmax, whose value is Mmax 



A is function only of r; it is given by: 
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If the E value is positive this A(r) function may have 
two minima. The inner minimum appears when the A value 
increases. When the total energy of the gas is negative we 
expect the system to be bound and therefore the solutions 
describing gas inflow into the black hole should not extend to 
infinity. Now the A{r) function goes asymptotically to zero 
at a maximum radius and may have a relative minimum, at 
the radial position ri, close to the black hole if the angular 
momentum is large enough. Figure 1 shows the A function 
for A = 1.67 and different E values (-0.00001, -0.00003, - 
0.00009, -0.00027); higher lines correspond to larger E values 
(smaller absolute values) . Figure 2 shows the A function for 
the same energy E = —0.00001 and different, increasing, 
angular momentum values (1.67, 1.67 x 1.01 , 1.67 x 1.02 , 
1.67 X 1.03); higher lines correspond to lower A values. Both 
figures are drawn with S = Q and r starting from r = 1.5 
to avoid the singularity at r = 1. For the same S value and 
with 7 = I we have the value w = 3.181980514. 

The radius ro at which A[r) goes to zero can be found 
analytically; it depends only on E and A, not on the 7 of the 
gas. 

Since the product of the Mach function f{M) and the 
A{r) function must be constant along the flow and since the 
Mach function has a single maximum, then the minimum of 
A{r) at ri must be coincident with the maximum of f{M). 
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Figure 1. A(r) function for different E values. Higher lines cor- 
respond to larger E values (smaller absolute values). 




Figure 2. A(r) function for different angular momentum values. 
Higher lines correspond to lower A values. 



Imposing f{M)A{r) = f{Mmax)A{r{) we obtain the 
solutions for M. Starting from the value f{Mmax)A{r\) we 
may obtain a new variety of solutions just solving for M the 
implicit equation f{M)A{r) = f{Mrr^ax)A{rl)^,. The factor ^ 
scales the entropy of the flow. Decreasing the value of ^ < 1 
we obtain families of solutions with the Mach number M 
shifted to larger values. 

For any supersonic solution Msup, by using the Hugo- 
niot relation we may calculate the post shock Mach value if 
a shock at a generic radius value r occurred: 



M, 



post 



up ^ 

7-1 



(8) 



A standing shock can occur in the solution at Vshock if 
the values Mpost{r shock) and Msubir shock) are equal: 



Msubifshock) — Mpost{T shock) 



(9) 



Among the solutions with different ^ values (i.e. with 
different entropy values) we will find the ones that can be 
connected by a shock jump with the subsonic branch. Figure 
3 shows the solutions for the Mach number with 7 = 4/3, 
A = 1.67 and E — —0.00001. The inner curve corresponds to 
5 = 1, the larger one to ^ = 0.5. It is clear that the solutions 




Figure 3. Analytical solution with shock. Mach number versus 
radius. 



can reach a maximum radius rm = 21.75: they do not extend 
to infinity. For this case the shock is at r shock = 6.85. 



3 NUMERICAL SIMULATIONS 

We report here the results of simulations performed us- 
ing a Smoothed Particles Hydrodynamic (SPH) code in 
cylindrical coordinates developed by the authors. The 
SPH formulae have been repo rted in previous publications 
JChakrabarti fc Moltenilll993^ . We stress here that a very 
accurate conservation of the angular momentum of the flow 
is an essential condition to obtain the shock solutions. For 
the numerical simulations we use dimensionless quantities. 
We adopt the speed of light as reference speed and the grav- 
itational radius of a Schwarzschild black hole Vg — 2GM* /c? 
as reference length. 

The comparison of simulation data with 2D analytical 
solutions requires some care. In the analytical treatment we 
assumed vertical equilibrium, but numerical simulations al- 
low the gas motion in the vertical direction. Therefore it 
is quite reasonable that the shocks, in the simulations, will 
not be exactly coincident with the theoretical predicted po- 
sitions. We add, as comment, that we made simulations also 
with a flctitious damping acting only on the speed. The 
vertically damped simulations show a shock position in bet- 
ter agreement with theory. 

We present here the data coming from the simulation 
of a 2D case with the parameters used to obtain Figure 
3. This case has the following parameters: E — —0.00001, 
A = 1.67. The shock is predicted to be at Vshock ~ 6.85. We 
chose to inject matter at ri„ = 20. With this parameters the 
radial speed and the sound speed at ri„ are Vi„ = —0.09606, 
din = 0.07790 respectively. The spatial resolution of the 
SPH particle is ft = 0.2. The steady state particle number 
is Np = 63000. 

Figure 4 shows the particle distribution in the up- 
per half of the domain. The shock of the simulation is at 

Tshock-SPH = 8.7. 

Taking into account our previous remarks we may say 
that the shock position of the 2D simulation is in good 
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Figure 4. 2D simulation: Particle positions. 
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Figure 5. ID simulation: Mach versus radial distance. The con- 
tinuous line is the theoretical one. 



agreement with the theory. We also report in Figure 5 
the Mach value versus the radial distance of a ID simu- 
lation. The physical parameters for the ID simulation are 
E — —0.003, A = 1.85, 7 — 5/3; the numerical parameters 
are ri„ = 75.26, = -0.04719, ai„ = 0.03931, h = 0.125. 
In this case, since there is no vertical motion, the agreement 
between theory and simulation must be exact, and effec- 
tively this occurs. For clarity, in the figure we plotted the 
theoretical Mach value plus a small shift(0.05). The shock is 
at the predicted position of Vshock = 12.8. The whole shape 
of the curve profile is very well reproduced. Note that the 
ID solution is obtained in the same way of th e 2D case, but 
the f unctions A{r) and f{M) are different. iMolteni et alJ 



4 CONCLUSIONS 

With our analysis we reinforce the view that the existence 
of permanent shocks around black holes is an important 
aspect of the physics of accretion flows. It has been sug- 



gested that the QPO oscillations of the X-ray luminosity of 
galactic black hole candidates m ay arise from shock oscil- 
lations dChakrabarti et al]|2004) .Jets of out-fiowing matter 
may arise als o from the high energetic po st shock gas driven 
by radiation iChattopadhvav et alJ2o'Ml . These shocks pro- 
duce around the central object a bulge of hot gas comptoniz- 
ing the low energy thermal emission from a cold keplerian 
disk. We remark that shocks produce the hot corona in a 
quite natural way, while the corona outcome or the ADAF 
solutions originating from keplerian disks have been show n 
by Molteni et al. to be very unlikely jMolteni et alj|200lll . 

We stress that these fiows are essentially axisymmetric. 
So this situation is different from the one in which the gas 
fiow, in a binary system, is captured from the wind origi- 
nated by the normal star. Steady shocks of this kind cannot 
be produced in this case since the flow has low angular mo- 
mentum, but has non-axisymmetric configuration. 

Obviously a crucial point is the question of the origin 
of such symmetric sub-keplerian flows. It seems to us that 
this problem is not critical for black holes in the centers of 
galaxies or AGN, since in this case the star environment is 
nearly spherically symm etric and with low global rotation 
JCohen fc Rvzhovl 119971) . In the case of black holes in bi- 
nary systems we suggest that accretion may occur not only 
from the inner lagrangian point, but also by a small amount 
of ambient gas, possessing low rotation, in the potential 
well o f the binary. Time dependent 3D numerical simula- 
tions iMolteni et alj|200ll) and recent observations support 
this hvpothesis. see lSniitj]_et_jJj 
suggested by CarrolJCano^D~ 



g002'). Furthermore it was 
L, 1990) that the Poynting- 



Roberston effect can be efficient also to drive matter from 
keplerian disks even in accretion on White Dwarfs. In the 
case of accre tion on weakly magne tized neutron stars it has 
been shown jMiUer fc Lamblll993l) that if the luminosity of 
the star is ~ 0.2Le { Le is the Eddington luminosity) a sub- 
stantial fraction of the accreting matter can lose its angular 
momentum. We also think that, for the binary black hole 
case, the disk self-illumination may produce similar effect 
decreasing the angular momentum of the outer layers of the 
disks. This process would be similar t o the accretion driven 
by ra diation suggested by Ballantyne jBallantvne fc Everetd 
|2005^ to occur in neutron star luminosity bursts. We add a 
trivial consideration: even if the inflow parameters do not 
correspond to the exact analytical ones, but are close to 
them, steady shocks may be still formed. 

If requested, the authors can offer the software to re- 
produce both the analytical and the numerical results. 
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